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We present results of a search for D^'-D" mixing and a measurement of _Rd, the ratio of doubly 
Cabibbo-suppressed decays to Cabibbo-favored decays, based on an analysis of Z?" K^-k~ decays 
in 57.1 fb~^ of data collected at or just below the T(4S) resonance with the BABAR detector at 
the PEP-II collider. Our results are compatible with no mixing and no CP violation. At the 95% 
confidence level, allowing for CP violation, we find the mixing parameters x''^ < 0.0022 and —0.056 < 
y' < 0.039, and the mixing rate Rm < 0.16%. In the limit of no mixing, Ru = (0.357±0.022 (stat.)± 
0.027 (syst.))% and the CP-violating asymmetry Ad = 0.095 ± 0.061 (stat.) ± 0.083 (syst.). 

PACS numbers: 13.25.Ft, 12.15.Ff, 11.30.Er 

Within the Standard Model the level of mixing ments fH . For this reason D^-D^ mixing is a good place 

is predicted to be below the sensitivity of current experi- to look for signals of new physics beyond the Standard 
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Model 0- Because new physics may not conserve CP, 
it is important to consider CP violation when measuring 
mixing. Observation of CP violation in D^-D^ mixing 
would be an unambiguous sign of new physics jQ] ■ 

Mixing can be characterized by the two parameters 
X = Am/r and y = Ar/2r, where Am — mi — m2 
(Ar = Fi — is the difference in mass (width) between 
the two mass eigenstates and F is the average width. 

The dominant two-body decay of the is the right- 
sign (RS) Cabibbo-favored (CF) decay D° ^ K-'it+. 
Evidence for mixing and CP violation, if present, will 
appear in the wrong-sign (WS) decay D'^ — s- K^i:^. 
Charge conjugates are implied unless otherwise stated. 
Two amplitudes contribute to the production of this fi- 
nal state: the tree-level amplitude for doubly Cabibbo- 
suppressed (DCS) decay of the D°, and an amplitude for 
mixing followed by CF decay of the 13° . Assuming that 
X, y <C 1 and CP is conserved, and with the convention 
AF ^ T{CP +1) - F(CP = -1), the time-dependent, 
WS decav rate Tws{t) for D° K^n~ can be approxi- 
mately y related to the RS decay rate TRs(i) by 




(1) 



In Eq. t is the proper time of the _D" decay mea- 
sured in units of the lifetime tjjo, T'Rs(t) oc e~*, 
i?D is the time-integrated rate of the direct DCS de- 
cay K^TT^ relative to the RS decay, and x' , y' 
are related to x, y hy x' — xcosSktt + ysinSKTv and 
y' = —xsihSk-tt + ycosSK-w, where Sktt is the rela- 
tive strong phase between the CF and DCS amplitudes. 
Physics beyond the Standard Model may include addi- 
tional phases that are not CP-conserving. Such terms 
can be absorbed into a phase ip, described below. The 
time-integrated WS decay rate is 



^ws — 



(2) 



Previous experiments have searched for mixing using 
wrong-sign hadronic |^ 1^ and semileptonic 7] 
decays, or have searched for width differences between 
CP = -Hi and CP = -1 states directly .S, JtllOl]- Since 
x' appears only quadratically in Eq. QJ, its sign can- 
not be determined in an analysis based on the WS decay 
alone. 

To allow for CP violation, we apply Eq. Q to D'^ and 



separately. We determine {R 
candidates and {R 



:+S,_x'+^ y'+} for 
' , y'~} for candidates. The 
separate and results can be combined to form the 
quantities 



WS' 



Pp - i?D 
— — T 



P+ 

A — M 



R 



M 



+ R 



(3) 



M 



where P^ = [x'^ -\- y'^ )/2. Ad and Ayi are related to 
CP violation in the DCS decay and mixing amplitudes. 



respectively. CP violation in the interference of DCS 
decay and mixing is parameterized by the phase Lp: 



1± A 



M 



itA 



M 



(x' cos(/3 ± sin(^), (4) 



1± A 



M 



(2/'cos(/3Ta;'sin(/3). (5) 

J- T 

An offset in of ±7r can be absorbed by a change in 
sign of both x' and y', effectively swapping the definition 
of the two physical D° states without any other observ- 
able consequence. To avoid this ambiguity, we use the 
convention that \if \ < 7r/2. 

We select a very clean sample of RS and WS decays and 
fit for signal and background components in a 57.1 fb~^ 
dataset collected with the BABAR detector at the 
PEP-II e+e^ storage ring. We extract the parameters de- 
scribing mixing and DCS amplitudes from the WS decay- 
time distribution. To avoid potential bias, we finalized 
our data selection criteria and the procedures for fitting 
and extracting the statistical limits without examining 
the mixing results. 

We select D° candidates from reconstructed 
D^ir'^ decays; this provides a clean sample of decays, 
and the charge of the pion (the 'tagging pion') identifies 
the production flavor of the neutral D. We retain each 
RS and WS Z?" candidate whose invariant mass niKiv is 
within 60 MeV/c^ of the mass. We require the mass 
difference Sm between the D*'^ and the D° candidate to 
be less than m^-|-25 MeV/c^. Only D*+ candidates with 
center-of-mass momenta above 2.6 GeV/c are retained, 
thereby rejecting D*~^ candidates from B decays. 

We determine the vertex by requiring that the 
D'^ decay tracks originate from a common point with 
a probability p(x^) > 1%, and then determine the D*^ 
vertex by extrapolating the D'^ flight path back to the 
beam-beam interaction region. This procedure benefits 
from the small vertical size (« 7 /xm) of the luminous 
region and the well-measured decay products. We 
constrain the trajectory of the tagging pion to originate 
from the D*~^ vertex, thus improving the measurement of 
Sm. We then calculate the proper time t of the decay 
from the dot product of the D° momentum vector and 
flight vector, defined by the D*+ and D° decay vertices 
in three dimensions. The typical resolution is 0.2 ps. 

We determine the mixing parameters by unbinned, ex- 
tended maximum-likelihood fits to the RS and WS sam- 
ples simultaneously. We perform four separate fit cases: 
first, a general fit allowing for possible CP violation, 
which treats WS and 13° candidates separately, fitting 



for {P+ 



wsi 



y'"*"} for candidates and {R 



WS' 



y'^} for candidates; second, a fit assuming CP con- 
servation, which does not differentiate between and 
candidates, fitting for {Pws, x''^, y'}; third, a fit as- 
suming no mixing, but allowing CP violation in the DCS 
amplitudes, fitting for {Pd, Ad}; and fourth, a fit for 
Pd, only, assuming CP conservation and no mixing. 
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For each fit case we assign each candidate to one of four 
categories based on its origin as I?° or D", and its decay 
as RS or WS. For each category we construct probabihty 
density functions (PDFs) that model signal and back- 
ground components. As independent input variables in 
the PDFs we use the candidate mass niK-w, the mass 
difference dm, and the D'^ proper time t with its error at- 

The fit is performed by simultaneously maximizing in- 
dividual extended likelihood functions, one for each can- 
didate category. Within each category, the likelihood is 
a sum of PDFs, one for each signal or background com- 
ponent, weighted by the number of events for that com- 
ponent. Each component's PDF factorizes into a portion 
describing the behavior of each independent variable con- 
voluted with a corresponding resolution function. The 
parameters describing the mass resolutions and shapes 
and the lifetime resolution are shared between PDFs. 
These parameters are determined primarily by the much 
larger RS sample. 

We characterize the WS background by three compo- 
nents: true Z)° decays that are combined with unas- 
sociated pions to form D*+ candidates; combinatorial 
background where one or both of the tracks in the D° 
candidate do not originate from decay; and back- 
ground where the kaon and the pion in a decay 
have both been misidentified, thus converting a RS de- 
cay into an apparent WS decay (double misidentifica- 
tion). Kaons (pions) are identified with an average effi- 
ciency of 84% (85%); the average misidentification rate is 
3% (2%). Correctly fitting the WS double misidentified 
background is particularly important due to the large size 
of the RS sample; its level as obtained from the fit agrees 
well with predictions based on our particle identification 
performance. 

We treat the normalization of WS candidates originat- 
ing as a D'^ or separately, thus yielding in total two 
signal and six background components in the WS part of 
the PDF. We assume CP conservation in the RS data; its 
PDF has one signal and three background components. 

We perform each fit in steps. Parameters correspond- 
ing to the mxTi and Sm distributions and the number of 
candidates in each category are determined first. In a 
second step, these parameters are fixed and a fit to the 
proper time distribution is performed. The shapes of the 
distributions in rriKTr and dm allow the fit to differentiate 
between the various signal and background components. 
Figure shows projections from the WS sample of the 
jTiKTr and Sm distributions overlaid on the fit result. 

We fit the RS decay-time distribution using a model 
that combines the RS signal decay-time distribution 
{TRs{t) in Eq. ^) and the expected decay-time distri- 
butions of each background component, convolving each 
with a common decay-time resolution model that uses 
the decay-time error for each candidate and a scaling fac- 
tor determined in the fit. For the WS signal component 
we use the same resolution model but with a lifetime 



distribution including the mixing parameters as given 
by Tws(i) in Eq. Q or its CP- violating counterparts. 
For the unassociated pion and double misidentification 
backgrounds we also use the TRs(i) lifetime distribution 
because they are true decays. The combinatorial 
background is assigned a zero-lifetime distribution and 
a signal-type resolution model based on studies of mass 
sidebands and Monte Carlo (MC) samples. 




FIG. 1: The distribution of the WS data for a) ttikit with 
144.5 < 5m < 146.5 MeV/c^, b) Sm with \mK7v - m^yol < 
20 MeV/c^ c) mx^ with 150 < Sm < 165 MeV/c^, and d) 
Sm with 25 < Imx-w — ?ti_do| < 60 MeV/c^. Data are shown 
as points with the contributions from the fit overlaid: signal 
(open), unassociated pion background (dark shaded), dou- 
ble misidentification background (black), and combinatorial 
background (light shaded). 



TABLE I: Fit parameter results determined by the full fit, 
with no constraint on x''^ in the mixing- allowed cases. For 
the no-mixing cases, R^l = R\^\ The +(-) signifies D° 



Fit case Parameter Fit result (/lO '^) 









0« 


D° + D° 


Mixing 
allowed 


r(±) 

-"ws 
:r'(±'' 


3.9 
-0.79 
17 


3.2 
-0.17 
12 


3.6 
-0.32 
13 


No mixing 


p(±) 

-"-WS 


3.9 


3.2 


3.6 



In Tabledwe summarize the central values returned by 
the fit for the four cases. In Fig.|21we show the decay-time 
distribution of the WS sample for the signal and a back- 
ground region. We select a signal (background) region 
to provide a sample with 73% signal (50% combinatorial 
background) candidates based on the reconstructed val- 
ues of m^Tr and Sm. The selected signal region contains 
64% of all signal events according to the fit. In total we 
observe about 120,000 RS (430 WS) signal decays. 

Our fit permits x'^ to take unphysical negative values. 
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FIG. 2: The proper time distribution for the WS candidates 
in a) the signal region (73% signal purity) and b) a back- 
ground region (50% combinatorial background). See Fig. 
for component definitions. 



The interpretation of non-physical results and error es- 
timates calculated from the log-likelihood surface (LLS) 
would require a Bayesian analysis where the choice of 
prior is not clear. In addition, an accurate error estimate 
from the LLS requires a LLS shape that is not strongly 
dependent on the outcome of the fit. These requirements 
are not satisfied here. Therefore, we use a frequentist 
approach, and construct 95% confidence-level (CL) con- 
tours in (x' , y') utilizing toy MC experiments. In each 
toy MC experiment we generate a WS dataset (the part 
sensitive to mixing) for a given (x' , y') with the same 
number of I?° and I?° events as observed in the data, 
but with a decay-time distribution appropriate for the 
chosen point. Fit parameters for the rriKTr and Sm dis- 
tributions and other parameters not sensitive to mixing 
are fixed at their fitted values from data. The at dis- 
tribution and background fractions from the data fit are 
used as well. We fit each toy MC dataset, obtaining 
values for the mixing parameters and the correspond- 
ing LLS. We construct contours such that for any point 
ttc = {x'1, y'c) on the contour 95% of the experiments gen- 
erated at that point will have a log-likelihood difference 
A In £(ac) = ln£inax"ln£(ac) less than the correspond- 
ing value Aln£data(ac) evaluated for the data. £max is 
the maximum likelihood obtained from a fit to either data 
or a toy MC sample. 

Where we assume CP conservation we apply this 
method to the combined D° and WS samples. The 
resulting contour is shown by the dotted line in Fig. 13 
The 95% CL limits for Ru and for Rm are obtained by 
finding their extreme values on the 95% CL contour. 

To consider CP violation, we divide the WS sample 
into candidates produced as a D° or as a and cal- 



culate separate contours for {x'^ ,y'^) and {x' 



_2 



,y' 



each corresponding to a CL of 1 — \/Q.05 — 77.6%. Each 
point on the D*^ contour is combined with each point on 
the D° contour using Eqs. ©-© to produce two poten- 
tial solutions of {x''] y'} for each relative sign of x''^ and 
x'~ . The outer envelope of these points is presented as 
the 95% CL contour in the (x'^ y') plane (see Fig.OJ. The 
peculiar shape of the contour arises from the two poten- 
tial solutions for each point on the Z?" and Z?" contours. 




O Physical (x , y') 

• Central (x^.y'j No CPV 

95% CL CPV allowed 

95% CL CPV allowed, slat only 

- - 95% CL CP conserved 

■ ■ ■ ■ 95% CLCP conserved, stat only 
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FIG. 3: 95% CL hmits in ix'^,y') with and without CP vi- 
olation (CPV) allowed. The sohd point represents the most 
hkely fit point assuming CP conservation and the open circle 
the same but allowing CP violation and forcing x'^ > 0. The 
dotted (dashed) line is the statistical (statistical and system- 
atic) contour for the case where no CP violation is allowed. 
The solid and dash-dotted lines are for the corresponding case 
where CP violation is allowed. 



This contour is more stringent than the CP-conserving 
case in some cases, which is allowed as the definition of 
coverage is slightly different. No central value for x'^ 
exists if either a;'+ or x'~ < 0. 

We summarize results for all four fit cases in Table ITTl 
We obtain limits on the individual mixing parameters by 
projecting the contours onto the corresponding coordi- 
nate axes. Since the no-mixing solution is well within the 
95% CL contour, we cannot place limits on Au and (p. 



TABLE II: A summary of our results including systematic 
errors. A central value is reported for the full fit with x'^ 
fixed at zero. The 95% CL limits are for the case where x' 
was not constrained during the fit. 



Fit case 


Parameter Central value 


95% CL interval 






{x"=0) (/10-3) 


(/lo-^) 


CP 

violation 
allowed 


Rd 

Ad 
a 

X 

y' 

Rm 


3.1 
1.2 


8.0 


2.3 < i?D < 5.2 
-2.8 < Ad < 4.9 

x'^ < 2.2 
-56 < y' < 39 

-Rm < 1.6 


No CP 
violation 


Rd 

/2 

X 

y' 

Rm 


3.1 


8.0 


2.4 <Rd < 4.9 
x'^ < 2.0 
-27 <y' <22 
Rm < 1.3 


No mixing 


Rd = 
Ad = 


(0.357 ± 0.022 (stat 
0.095 ±0.061 (stat.) 


) ±0.027 (syst.))% 
±0.083 (syst.) 


No CP viol, 
or mixing 


Rd ~ 


(0.359±0.020 (stat.) 


±0.027 (syst.))% 
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To estimate systematic uncertainties we evaluate the 
contributions from uncertainties in the parametrization 
of the PDFs, detector effects, and event selection criteria. 
The small systematic effects of fixing the rriK-K and 5m 
parameters and the number of events in each category 
in the final fit is evaluated by varying these parameters 
within statistical uncertainties while accounting for sta- 
tistical correlations. 

For detector effects such as alignment errors or charge 
asymmetries we measure their effect on the RS sample. 
Under the assumption that RS decay is exponential and 
has no direct CP violation, this method is very sensitive. 
The systematic error due to the size of the MC sample is 
insignificant since all distributions are obtained from the 
data. 

Each systematic check yields a small shift in the fitted 
mixing parameters. We use MC experiments to deter- 
mine the significance of each shift using the same method 
employed as for the 95% CL statistical contour. We scale 
the statistical contour with respect to the central fitted 
point by the factor -y/ 1 -I- ^ m| , where is the rela- 
tive significance of each systematic check. For the gen- 
eral case we carry out this procedure for the and 
contours separately before combination. In all fits the 
largest effect for x'"^ and y' is the D*'^ momentum se- 
lection cut, with mf = 0.24; all others are at least three 
times smaller. For i?D the largest effect is the decay-time 
range. We show contours including systematic errors in 
Fig. O as a dashed line in the CP conserving case and as 
a dash-dotted line in the general case. 

In summary, we have set improved limits on D^-D'^ 
mixing and on CP violation in WS decays of neutral D 
mesons. Our results are compatible with previous mea- 
surements ^ 1^ and with no mixing and no CP viola- 
tion, which agrees with Standard Model predictions. 
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